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has  befcn  notified  to  the  contrary  by  the  sponsor.  If  this  protective  statement  is  j 
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MASS  EARTH  MOVEMENT  AND 
TREATMENT  FEASIBILITY  STUDY 


1 


by 

Emmericus  C.W.A.  Geuze 

Professor  of  Soil  Mechanics  and  Foundation  Engineering 
Rensselaer  Polytechnic  Institute 

SYNOPSIS 

The  basic  concept  of  desaggregation  is  to  develop  a state  of  failure  in 
soil  media  by  a centre  of  pressure  located  at  a convenient  depth  below  the  sur- 
face of  the  soil  mass. 

The  mechanism  of  failure  development  in  the  medium  located  between 
the  centre  of  pressure  and  the  soil  surface  has  been  based  upon  the  results 
obtained  by  application  of  the  theory  of  plasticity  of  granular  and  cohesive 
materials. 

The  failure  properties  of  these  two  soil  groups  have  been  expressed  by 
the  parameters  of  internal  friction  and  of  cohesion. 

The  basic  problem  is  presented  in  an  analytical  form  which  yields  approxi- 
mate solutions  for  the  pressures  requited  to  cause  failure  in  the  medium  at  diff- 
erent depths  of  the  centre  of  pressure  and  with  different  diameters  of  the  pressure 
cavity. 

- The  consequences  of  the  approximations  are  discussed  and  an  altogether 
new  theory  is  developed.  The  results  of  an  application  of  this  theory  to  the 
failure  of  a granular  material  are  shown  and  a comparison  with  results  obtained 
by  the  approximate  analytical  solution  is  made. 

Liquid  loading  characteristics  are  discussed  in  relation  to  the  compressi- 
bility and  permeability  of  the  soil  materials. 

A model  experiment  has  been  used  to  demonstrate  the  development  of  the 
failure  area  and  to  provide  pressure  data.  49398 


2 


Analysis  of  the  basic  problem 

The  pressure  cavity  is  supposed  to  have  the  shape  of  a sphere  with  radius 
ro*  centre  located  at  the  depth  D below  the  horizontal  boundary  of  the  soil 
medium.  Pressure  is  applied  to  the  wall  of  the  cavity  through  the  intermediary  of 
a thin  rubber  membrane.  The  cavity  is  completely  filled  with  liquid  (water),  which 
is  held  under  pressure  (Fig.  1). 

A certain  pressure  is  required  to  maintain  the  initial  diameter  of  the  cavity. 
An  increase  of  pressure  is  accompanied  by  outward  radial  compression  of  the  soil 
material  around  the  cavity.  The  compressive  strains  decrease  rather  rapidly  with 
increasing  distance  from  the  centre  of  the  cavity;  approximately  proportional  to  the 
one-third  power  of  that  distance. 

In  a dry,  or  moist,  granular  material  compressions  and  the  subsequent  dis- 
placement of  the  cavity  wall  occur  instantaneously.  Compaction  in  the  immediate 
vicinity  of  the  cavity  is  accompanied  by  the  formation  of  a dense  shell  of  material. 
When  the  material  is  saturated  with  water,  the  rate  of  compression  depends  on  the 
compressibility  and  permeability  characteristics  of  the  material  and  on  the  dura- 
tion of  the  imposed  loading  increment. 

Though  the  concept  "rate  of  consolidation"  usually  is  not  applied  to  gran- 
ular materials,  the  process  of  radial  compression  in  the  saturated  medium  described 
previously  does  not  differ  essentially  from  that  in  a saturated  cohesive  material. 

A rapid  build-up  of  the  cavity  pressure  holds  the  depth  of  penetration  of 
the  consolidation  effect  to  a minimum.  Failure  of  the  material  occurs  when  a cer- 
tain critical  value  of  the  cavity  pressure  is  attained.  The  failure  wave  propagates 
through  the  material  in  upward  direction  as  lopg  as  the  pressure  is  maintained  at 
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the  critical  value  until  it  reaches  the  upper  boundary  of  the  medium.  When  the 

'V 

cavity  pressure  is  kept  at  the  same  value  beyond  this  point,  the  body  of  soil 
within  the  failure  area  continues  to  move  in  upward  direction. 

When  the  impervious  membrane  is  left  out  water  will  penetrate  into  the 
medium  as  a result  of  the  increase  of  cavity  pressure.  In  a dry  medium  the  pene- 
tration of  water,  along  with  the  displacement  of  air  in  the  voids,  occurs  as  a 
capillary  saturation.  The  rate  of  penetration  depends  on  the  compressibility, 
the  permeability  and  the  capillary  suction  force  of  the  soil . In  a saturated 
medium  the  rate  of  penetration  depends  on  compressibility  and  permeability  only. 

The  mechanism  of  this  penetration  process  has  identical  consequences 
for  the  transfer  of  the  liquid  load  to  the  soil  medium  as  when  an  impervious 
membrane  would  separate  the  soil  from  the  liquid,  i.e.  , the  load  acts  upon  the 
medium  through  the  intermediary  of  a spherical  soil  shell.  This  shell  is  com- 
posed of  more  or  less  compacted  soil.  Its  thickness  depends  on  the  rate  of 

/ 

building  up  of  the  cavity  pressur  e and  on  the  soil  characteristics  as  mentioned  1 

earlier. 

Chapter  I 

The  approximate  analytical  solution  of  the  basic  problem 
An  approximate  analytical  solution  has  been  given  by  NADAI  [l]  for  the 
problem  of  a small  spherical  cavity  held  under  static  pressure,  situated  deep 

\ 

under  the  surface  of  granular  material. 

NADAI  solved  the  problem  by  making  two  simplifying  assumptions: 

a.  The  state  of  stress  in  the  vicinity  of  the  axis  OA  (Fig.  2)  is  taken 
approximately  equal  to  that  in  a spherical  body  of  the  same  material 
around  O,  tangent  to  the  surface  at  A. 
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b.  The  gravity  field  for  the  material,  which  has  the  direction  AO  at  every 
point  in  the  medium,  is  replaced  by  a field  of  body  forces  converging 
radially  toward  the  centre  O and  of  the  same  magnitude  as  the  uni- 
form force  of  gravity.  The  approximate  nature  of  Nadai's  solution  is 
a consequence  of  these  assumptions. 

The  equation  of  equilibrium  of  an  element  of  the  medium  in  the  shape  of  a 

truncated  circular  cone  depends  on  the  radius  (Fig.  3)  only,  because  of  the  polar 

symmetry  of  the  stress  field. 

Equilibrium  of  the  element  in  a radial  direction  requires  that: 

X do . x do  _ (G**  C*cx.)’(l+  dx.)do.(x+  dx. ) de  + 2.C^.  t.do+{i+oi\)c/oj  dr. 

-V.d  jzde.tde  + (x+dx)do  .{x+dx.)d<ijdx.  = C 

Dividing  through  by  't‘(do)X and  neglecting  ctxjx.  with  respect  to  unity: 

+r.x  ~ o (l) 

This  equation  can  be  solved  under  the  assumption  that  failure  occurs  simultaneous- 
ly along  sliplines  located  in  radial  planes  through  the  axis  OZ  and  in  tangential 
planes  through  O at  right  angles  to  the  radial  planes.  The  major  principal  stress 
0^  acts  in  radial  direction,  both  minor  principal  stresses  C£T  in  tangential 
directions. 

A.  Granular  material 

We  introduce  the  Mohr-Coulomb  condition  of  failure  in  the  equation  of 
equilibrium: 

CTx  ' ..C~ 

(2) 

Putting:  Vs'*  - a.  (a  parameter  of  the  solution) 

Eq.  (1)  takes  the  form: 
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The  solution  of  this  equation  proceeds  by  introducing  a new  variable: 

X = £L 

X 


hence 


is. 

X.  aCt» 


(4) 


and  Eq.  (3)  takes  the  form; 


and  as: 


Z^X  + O 

after  rearranging  terms  and  dividing  through  by  r: 

cL  X _ 

((aL+i)~)C  +>  * 

dX=fihrc>  d[^'i>c+r] 

/ a 

' [ccL+‘)y-+Y]  * 

Assuming  failure  to  occur  at  a cavity  pressure  pQ,  we  have  the  boundary  condition: 

hence:  /.*»*■ 

/ f d{?a.+i)X-  ^ 

X.+-irJ  J 'L 

-V*  u?t 

yielding  the  equation  , 

It, 

fr**')  X.W 

from  which  we  obtain:  {*■+>) 

ffa+o* +Y]  ~ y7^rj 

and  after  replacing  XbyCJT/i  and  Xftby  />»/*•  : 

ffl-fw-Aftr-Zr, 

Multiplying  through  by  r yields: 


*•  A jS* 


(5) 


The  validity  of  this  equation  is  limited  by  the  condition  that  the  state  failure  can 
only  exist  when  the  major  principal  stress  <T  t Q , because  the  Mohr-Coulomb 
failure  condition  is  expressed  by: 


according  to  Eq . (2) . 


Hence  the  external  boundary  of  the  spherical  medium,  when  free  of 
stresses,  represents  the  limit  of  the  medium  in  failure. 

This  condition  provides  the  other  boundary  condition  required  for  determining 
the  radius  of  the  spherical  body  of  failure. 

It  follows  that  in  order  to  find  the  magnitude  of  this  radius,  we  put  = O 
for  a radius  r^  in  Eq.  (3). 


(7) 
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If  the  quantity  is  given  besides  the  quantities  ex.  , Y and  r^ , the  cavity 
pressure  pQ  can  be  computed  from; 

r (<*-+•)  1 

-*•  ■ &-1®  - j 

Fig.  4 and  5 show  diagrams  of  Pq-^  relationships  for  rQ  = I,  5,  10  and 
.5  err.,  respectively.  The  volume  weight  of  the  material  was  taken  as  1 gram/cm'3 
and  the  angle  of  internal  friction  of  the  granular  material  y?  = 30°.  Hence  the 
coefficient  a = 4/3. 

It  should  be  noted  that  the  radius  of  the  cavity  has  a considerable  effect 
on  the  magnitude  of  pQ  since  it  appears  in  the  denominator  of  the  ratio  r^/rQ.  Gran- 
ular materials  are  characterized  by  an  angle  of  internal  friction  y?>  o.  Hence  the 
factor  (a  + 1)  has  values  larger  than  unity. 

Model  rules 

Eq.  (8)  is  of  significant  value  in  establishing  the  design  of  model  experi- 
ments as  will  be  shown  later. 

Since  both  the  volume  weight  X and  the  factor  (a  + 1)  can  be  adjusted  in 
the  experiment  according  to  prototype  conditions,  we  can  establish  the  following 


model  relationships 


where  capitals  indicate  the  relationship  for  the  prototype  and  the  small  characters 


(*-■+•)  - 

~p  _ r:7?c.  //^  ) - / 

° ~ A rJ  ! 

A ft)  ''I 


indicate  the  relationship  for  experimental  conditions. 


Setting: 


5.  - 

7?„ 
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we  obtain  the  ratio: 


or:  Z>  — ^ >b 

^ “ V ° (10) 

Eqs.  (9)  and  (10)  establish  the  model  rules  for  experiments  on  a small  scale, 
which  can  be  used  to  predict  prototype  behavior. 

Rate  of  dissipation  of  the  cavity  pressure  within  the  medium. 

The  rate  of  dissipation  of  the  cavity  pressure  is  expressed  by  Eq.  (5). 
When  dividing  through  by  the  following  equation  results: 


py  * M ® 


(<*~ +•) 


When  the  quantities  Y,  a,  rQ  and  pQ  are  given  constants  Eq.  (11)  can  be 
written  in  a simplified  form: 


(&■)*(*+ 


The  rate  of  dissipation  depends  mostly  on  the  magnitude  of  the  power  1/B,  i.e. , 
on  the  angle  of  internal  friction  <f>  . As  a demonstration  of  the  rate  of  dissipation 
we  take  V = 1 gram/cm^  and  - 19°26'. 

THe,  (5.)  - (A+±){±f-i 

i.e. , the  ratio  of  is  approximately  proportional  to  the  square  of  the  ratio 

V r. 

Fig.  6 shows  a diagram  of  the  <Tr- r relationships  when  2T=1,  rQ  = 10  cm, 
a=4/3  and  pQ  « 103,  104  and  5.104  grams/cm2. 

The  sllpplane  field 

The  slip  planes  intersect  the  principal  stress  trajectories  at  angles: 

fi,  * = f 

A=  f - 
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The  trajectories  of  major  principal  stress  are  a family  of  straight  lines 
radiating  from  the  centre  of  the  cavity.  The  trajectories  of  minor  principal  stress 
are  a family  of  concentric  spherical  surfaces. 

The  family  of  slip  planes  intersecting  the  family  of  radii  at  an  angle  of 
(y/+S)  are  called  Sj-planes,  those  at  an  angle  of  (y^-S ) are  called  S2-planes. 

The  s^-planes  and  s2-planes  are  conjugate  planes. 

Each  set  of  conjugate  planes  forms  a body  of  rotational  symmetry  with  re- 
spect to  a radius.  The  intersection  of  this  body  with  a plane  passing  through  the 
radius  is  a set  of  conjugate  slip  lines,  which  are  also  called  s^-  and  splines. 

The  derivation  of  the  Sj-  and  s2~slip  lines  follows  from  the  basic  property 
of  the  direction  of  slip  with  respect  to  the  direction  of  major  principal  stress  in 
the  medium . 

At  point  P in  the  diagram  of  Fig.  7 the  tangent  to  the  Sj-slip  line  has  the 
direction  (y'+S’)  and  to  the  s2-slip  line  the  direction  ty-S)  with  respect  to  the 
Z-axis. 


The  tangent  to  the  s^-slip  line  satisfies: 

= cotr,  Z 

7?.  bf 


Tailing  a point  RQ,  - y$,  at  the  boundary  of  the  cavity  (Fig.  8)  as  the  origin 
of  the  slip  line:  7? -7? 

f L5L  - cotn  r/V 

J *.  j 


7?,  R. 


V—J* 


Hence: 


and: 


Iff  JL  m Cotr* 

7?  s *•*/>  otr*  S".  +■  '■/'•) J 


(12) 


1 


i 

( 
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The  tangent  to  the  s2-sllp  line  at  P satisfies: 

^ - — cair-n  S' 

TK.'Zy 

because  the  sign  of  ^y>is  negative. 

Taking  a point  RQ,  + <f0  as  the  origin  of  the  s2-slip  line: 

^ — _ cotn  S’ j 

Hence:  _ _ r „ . . J (13) 

7Z  =R0.  <-*■/*[-  c°r”  *(+-%)] 

At  the  point  of  intersection  of  the  two  slip  lines  the  radii  R^  from  (12)  and 
(13)  are  the  same;  hence 

7?o.  /.j/*  /£.  **/>/- 

giving: 

y*c  « -y?+yi 

and:  , / 

and  by  substitution  in  either  Eq.  (12)  or  Eq.  (13),  we  obtain: 


= cjc/*fco+ri  r.  y>y 


(14) 


where: 

R^  * radius  of  the  failure  area. 

From  Eq.  (14)  follows  that  the  ratio  between  the  radius  of  failure  Rj  and  of 
the  cavity  Rq  increases  with  the  power  ro^rr  £”.  t/'  ♦ 

The  area  of  the  surface  of  a spherical  segment  with  radius  Rq  and  with  a 
centre  angle  zy£ represents  the  loaded  area. 

The  area  of  this  surface  is: 

SLTT  #*/',-  eo*  %) 

The  pressure  acting  on  this  area  can  be  obtained  from  Eq*  (8)  for  given 
magnitudes  of  Rc,  Rj  and  for  S~  * 30°. 

The  diameter  of  the  failure  area  is  significant,  because  it  establishes  the  degree 
of  approximation  involved  by  the  assumption  of  a radial  body  force. 


11 


From  Fig.  8 it  can  be  seen  that  the  strongest  deviation  of  the  radial  direc- 
tion from  the  axis  of  symmetry  (i.e. , the  direction  of  the  gravity  force)  occurs  at 
a point  where  the  tangent  to  either  the  s^  - or  s -slip  line  has  the  direction  of  the 


Z-axis. 


At  this  point: 


y/  - S"  = o 


“ / + 


5-  = o 


Hence: 


From  Eq.  (12)  follows  by  substitution; 

7?^  - 7?0.  ZKp j^cotn 

where;  R = the  radius  from  the  origin  to  the  tangent  point  on  the  slip  line. 

These  tangent  points  are  located  on  a circle  with  radius  rt,  where: 

r = ~R.  . s*  n S~ 

* * 


giving: 


in  S’.  Licp £cotn  y/jj 


and  using  Eq.  (14)  for  the  expression  of  Rj: 

= 7*,.  **/*[-  cotn  S'. 
we  obtain  by  substitution:  T5 

Z _ *inS.  mp [S'.  Co tnS] 

and:  r 7 

r.  CK.pj^S’.  CctnTj  (15) 

Fora  given  granular  material,  S’  is  a constant.  The  maximum  diameter  of  the 

* 

failure  area  is  therefore  proportional  to  the  depth  of  the  cavity.  The  radial  direc- 
tion of  the  body  forces  at  the  level  of  maximum  diameter  is  independent  of  the  depth 


of  the  cavity,  because: 


^ = Sin  r 


NADAI’s  approximate  solution  is  therefore  subject  to  a serious  deviation 
from  actual  prototype  conditions , where  the  direction  of  the  body  force  at  the 
level  of  maximum  diameter  may  differ  from  the  direction  of  gravity  as  much  as 
by  an  angle  £”=  (4  5°  - . This  deviation  decreases  with  increasing  magni- 

tude of  if.  i.e. , of  the  angle  of  internal  friction  of  the  material. 


i 
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5.  Cohesive  material 

In  a cohesive  medium  with  no  internal  friction,  y?  = o,  which  leaves  us  with 
only  one  parameter  for  the  failure  condition  of  the  material,  i.e. , the  cohesion  c. 
The  minor  principal  stress  ^ then  obtains  the  value: 

<£  * 

Substitution  of  this  expression  in  the  equation  of  radial  equilibrium  (1)  yields: 


VC  = o 

,,  - cit. 


d.C'  = -Y.dx-Uc  — 

z 


With  (T  9 h at  z,  - z : 

/C  /J.CJ  ✓ 

J cLc^  - I cL*.  - VC  j 


tenoe:  <T  , /,  _ r 

At  r = r , G~v-  o,  yielding: 

A = 

In  dimensionless  quantities: 
The  slip  line  field 


(17) 

o,  yielding: 

A = >7*,-V+VC./r.^) 

(18) 

quantities: 


The  equations  of  the  slip  lines  are  directly  obtained  from  Eqs.  (12)  and 
(13),  by  putting  (P  = o,  giving: 

hence:  co  tn  a * / 

and  the  Eqs.  change  to: 

and  / , t 

for  s - lines  and  s -lines  respectively. 

1 2 

The  point  of  intersection  follows  from  the  equality  of  Rj  at^ j; 

+, 


(21) 


and  by  substitution  we  obtain: 

Other  relationships  as  previously  obtained; 

Vt*  ° 

\ (*u  «-  vr*j 

Q,  d'  V7^ 

The  same  objections  as  previously  raised  apply  to  a cohesive  material.  The  direc- 
tion of  the  radial  body  force  deviates  by  as  much  as  45°  from  the  actual  direction 
of  gravity  at  the  level  of  maximum  diameter  of  the  failure  area. 

Conclusions 

The  analysis  of  an  approximate  solution  following  Nadai's  approach 
yields  results  which  may  show  reasonably  correct  values  in  the  vicinity  of  the 
axis  of  symmetry  with  a direction  parallel  to  the  force  of  gravity. 

The  slip  line  field  extends  to  a considerable  distance  from  the  axis  of 
symmetry  for  both  granular  and  cohesive  materials.  At  this  distance  the  direction 
of  a body  force  in  the  radial  direction,  as  assumed  in  the  approximate  solution, 
shows  a considerable  deviation  from  the  direction  of  the  force  of  gravity. 

A solution  as  suggested  by  Nadai  cannot  be  taken  as  a basis  for  the  con- 
venient computation  of  the  quantities  required  before  it  proves  its  validity  in 
comparison  with  the  results  of  an  exact  solution  of  the  problem. 


/r 


Chapter  II 

A.  The  failure  of  the  soli  by  cavity  pressure,  when  the  medium  Is  subjected 
to  a uniform  held  of  gravity  forces. 

The  basic  assumptions  are  identical  with  those  made  in  the  derivations  of  the 
preceding  analysis. 

The  stress  field  is  symmetrical  with  respect  to  the  direction  of  the  gravity 
force.  The  vertical  axis  is  taken  opposite  to  that  direction. 

Using  cylindrical  coordinates,  the  equilibrium  equations  can  be  derived 
from  the  equilibrium  of  a cylindrical  element  (Fig.  9).  Since  the  stress  field  is 
symmetrical  to  OZ,  radial  planes  are  subjected  to  a normal  stress  only.  Shear 
stresses  act  on  the  planes  r = constant  and  z = constant. 

We  have  for  equilibrium  of  forces  in  Z-direction: 

fc  -ft +15 <<*)]■  cU) c/e.tU 

7 _ yr  otT.'jcl&.ctl..  djt  —O 

* Ui.de.dh  J^L.dt.cLx.cia.db 


± /JLz  *.oix).do.ctzd'h  —T.ctx.do.cCk  - 


Dividing  through  by  r.da.dt.  di  gives: 

. - x 0+  $)  - *(•*  r)  - ° 

Since  we  may  neglect  cL'c/x.  with  respect  to  unity: 

* TZ 


(22) 


Equilibrium  of  forces  in  r-direction: 


C T.icl&.d*  Ui^i+dxjc/o.dt  + s.n  ^cta.  c/i.c/t.  + 

j-ft.'z.i-dij.d&.d'c  =.  o 


and  with; 


0>  / -X7  J.  c/q  a J.  d/o 
X X 


X.  c/e.  c/ 1.  c/Jk  — <.V«9.  «a^t.  c/*  _ c/r.  a/o.o/t.djt  -j-  C^.o/e.  c/t.  ci 

J'C 


+ Mc/iJ  c/e.  dx.c/yk 


Dividing  through  by  't.c/e.dt  dl  yields: 


and; 


-If-f -?f • £ + f - w/'**- r)  “ ° 

f - f+ If  ° 


Neglecting  ^2:  with  respect  to  unity: 

^ ^ /<z-<Z)  = o 

y* +’  x y* 

Eqs.  (22)  and  (23)  can  be  simplified  by  the  following  substitutions: 


'L.  /t.r) 


yielding; 


=V'^  ■*  " rr 


, r 


- _ ytx. 


(23) 


(24) 


h-^)*h^T}  = * 


(25) 


a 


The  Mohr-Coulomb  concept  of  failure  requires  tangency  of  the  stress  circles 
to  the  failure  envelope  (Fig.  10). 

The  relation  between  the  stresses  then  is: 


Introduction  of  the  failure  condition  (26)  into  the  differential  equations  of 

equilibrium  are  simplified  by  following  Nadai's  approach  [0  ■ 

The  stress  , (Th  andTTare  expressed  in  terms  of  the  variables  y and  ; 

^ s angle  of  the  r-axis  with  the  positive  direction  of  the  major  principal 
stress 

1/2  (<7+C£)  = mean  value  of  the  radial  and  axial  normal  stresses. 

From  Fig.  10  follows  that: 

(£«  <77,  (/+  eomf) 

5.  (27) 

r - fji,  ( Jim  </■  OH.  M li/J 

According  to  the  Mohr-Coulomb  concept  the  intermediate  principal  stress 

does  not  affect  the  ratio  between  the  major  and  the  minor  principal  stress  at 

failure.  Since  we  assumed  a failure  field  of  rotational  symmetry  with  respect  to 

the  z-axis,  the  Intermediate  principal  stress  is  the  tangential  normal  stress  acting 
0 

on  radial  planes  through  OZ.  The  tangential  normal  stress  can  therefore  obtain  any 
value  between  the  principal  stresses  (7^ and  Vj  without  infringing  the  failure  con- 
ditions. Since  the  mean  value  of  CJ and  CJ’equals  the  mean  value  of  and  7^  : 


/ 8 


7«e  magnitudes  of  CJ"  and  Cj  can  be  obtained  from  the  first  and  second  Eq.  of  (27) 
by  taking  y-  = o: 

cr=  (j  -*■  -s<  *»  •/’) 

Sj  = l'~ 

The  tangential  stress  therefore  satisfies: 

^7  > £ * 

and  we  can  write: 

C~u  ( t + *)-C7r,  (28) 

where:  — / < °<  <*♦* 1 

Elimination  of  q,  C"  and  Z from  Eqs.  (24)  and  (25)  with  the  help  of  Eqs.  (27)  proceeds 
as  follows: 

Z.CJ= 

l.T  = '5lM  *’/' 

from  which  we  obtain  the  following  derivatives: 

~L  (XT)  =.  sinf.  sir,Xsj,.']L-(i.ty+  *■/*■  ■<** 

^ £ 

= (/-  J.»y. 

L/i.cg 

» r 

Introduction  of  these  expressions  in  Eqs.  (*v)  and  (M)  yields. 

0 

«„y.  «i„  if . *-/».«.)  4-  *‘"V  M 

+ (>-  s‘" *+  '2*  <29) 


- 

t X.[x  C^J.  soMf-i-"-'! 

■f  S<  n y?  Si-n  2^. . ^--  £l.  s<«y).  <<41  i x.^4  . ^0) 

These  equations  are  simplified  further  by  introduction  of  the  slip  lines  as  a 
reference  system  instead  of  the  t- A coordinates. 

When  the  coordinates  of  the  diagram  of  stresses  are  oriented  with  respect 
;o  the  physical  plane  in  such  a manner  that  the  positive  direction  of  thetf"^  axis 
coincides  with  the  positive  direction  of  the  r-axis  (Fig.  10) , the  direction  of  the 
major  principal  stress  at  a point  A in  the  physical  plane  is  identical  with  the 
vector  Qq P in  the  stress  diagram.  The  slip  line  directions  follow  from  the  loca- 
tion of  the  tangent  points  and  , representing  the  shearing  stresses  at  the 

slip  planes  in  A.  The  lines  C^P  and  Q2P  represent  the  shearing  directions  with 
respect  to  the  direction  of  major  principal  stress  QqP,  hence: 

angleQ  PQ  =(45°-  ^) 
i o 

angleQ2PQo  = (45°  -£) 

The  angle  enclosed  by  Q^P  and  the  positive  -axis  therefore  is: 

ft  3 f -rj 

and  by  Q2P  and  the  positive  cT-axis: 

fa  - * - lvr‘-  t) 

The  positive  directions  of  the  slip  lines  at  point  A in  the  physical  plane  are 
indicated  by  arrows.  The  slip  line  direction  of  will  be  called  the  s^-direction; 
the  direction  of^  the  -direction. 

) 


lo 


The  partial  derivatives  with  respect  to  these  directions  can  be  built  up 
according  to  the  chain  rule: 


where: 


3 

■Di 

•3  a 

5? 

= 

~ i)  t 

/ 

' ^ 

_ ■* 

. 3 

■>2- 

"«) 

" ;>* 

‘5Z 

*■  ^y3/ 

0 

J 

^lu- 

ll 

's'r73/ 

* *"/*» 

> 

K 

. Si  H 

Substitution  of  these  expressions  in  the  first  two  equations  yields; 


L.  - 


(31) 


In  order  to  obtain  Eqs.  (29)  and  (30)  in  a form  identical  to  the  Eqs.  (31) 
we  perform  the  following  operations: 

Multiply  Eq.  (29)  by  and  Eq.  (30)  by  («s<. ryj.)  and  add: 

|L  (/  + ^ 

■* 1 ,w  + ^ */**] 

+ + **»*+'  *°Y*J  * (32) 

_C£..<S/*A  - 


A/ 


Multiply  Eq.  (29)  by  (-'^)  and  Eq.  (30)  by  (▼■*<»/ \)  and  add: 

i__  ( 1 .Zlr.) / -v  y>.  S*  y>.  s,y».  dosfS^J 

_ l in  lif,.  Sinfa  4.  <*Oi  Ay-  .dm 

Sin  - (/-  st-nf.  C.as  iyj  ^VV/ 

4.  l/i.fc).  s*»if.'hk^eci*.f.si*(Zl  - stni<j..eoz/ij  = 

C ^sinfa +.Y:‘i.  aos/3,  (33) 

The  expressions  between  brackets  can  be  largely  simplified  by  the  use  of: 

/V  /+/v«-*-w0 
/V  4 

i3,*/3-.  * *y 

A -A  •Ir-'f) 

A * */  - A 
A * if-/3, 


giving 


-/'+ 

>7  y>. 

'S'  *»/**  + 

Si  n y>. 

J/rt  iy,  - 

d d Zt  ^Ze 

c 4.  <*y>  j 

t' 

v/»  Xy. 

- 

£«j  a y - 

Coi  /A 

- — 

£ o ${3t 

+/>- 

Stny.  di>i  tyj- 

do*  fa  - 

■y/ »» y • 

. ,y<  *j  /.if, . 

■soy3* 

=.  4.  <?0S 

V" 

Si-nfa 

4-  <s.»  *y- 

COS/J-  - 

o£/3a: 

c & i *y • 

a iy,/^x 

— +- 

Sm  y£f, 

From  Eq.  (33) 

in  terms 

-=•  f 

./  * # /2 

/4>ixy-J- 

JH*\fa  - 

y>. 

»i  a.y « 

■ lCifa 

+ *<™*y.  *<»/*,  + e**x+ •***(>,  *■  ■*•  **‘A 

-//-  ftp, + *»'/-**'**+•*"' ft*  - + £*4f‘ 

4.  tf.  nfa  - Sinlf.im  A » - Sir,A 

Substitution  of  these  expressions  in  Eq.  (32)  yields: 


_/c£. «/*y£  «*y\)  (34) 


% 


andinEq.  (33): 

la  S yjt, 0. 5/3,  J—/r.  -r  0‘-/3.  • _jL  /*  ■ Cj7 - 2/r^J.  4.  j < ry^  : 

+ -r  yn.  <?<> -,//^//)  (35) 

By  comparing  the  expressions  between  brackets  with  those  on  the  right 
hand  side  of  Eqs.  (31)  it  can  readily  be  seen  which  substitutions  can  be  made  in 
Eqs.  (34)  and  (35) . 

These  substitutions  yield: 

eosy.  » 4 toup.)  (37) 

To  obtain  these  equations  in  a form,  which  shows  that  the  s^-lines  and 

s -lines  are  the  characteristics  of  the  solution  we  introduce  the  new  variable: 

2 

X - . Ir, /z.C-m)  l38> 


giving: 


_ coif'  _J_ 2_/?.<r) 

IS,  “ *>'3/ 


~>JL  = ?aiJL . <-  . 2- 

/?•<£.)  ‘si 

and  fr.C^Js  £*/>  ( L f- 

Substitution  into  Eqs.  (36)  and  (37)  yields: 


5 /x  + u.)* (39) 

r/ 

i /V  . HO) 

~rJ~  3>‘"r 


JL  3 


If  at  a point  A the  values  X.  and  ^ of  the  variablesXand  y'are  known, 

Zq.  (39)  provides  us  with  information  about  the  increase  of  the  sum  ()d+y)  along 
the  s - slip  line.  The  individual  values  of Xandy  at  a point  C on  this  slip  line , 
r.owever,  cannot  be  obtained  from  this  information  only.  It  is  therefore  necessary 
that  the  values  of  X^and  are  known  at  a point  B located  on  the  s2-slip  line 
intersecting  the  s^-slip  line  at  C.  In  that  case  the  increase  of  the  difference 
()C-y)  provides  us  with  the  information  about  the  increase  along  the  s2~slip  line 
when  going  from  B to  C,  using  Eq.  (40) . 

By  combining  the  information  about  the  increases  of  (X  + y')  and  (Xry) 
along  the  s^-slip  line  and  the  s2~slip  line  respectively,  X^and  y^can  be  computed 
individually  as  follows: 

(*c+ {*»  + *>)  + 4 

C^~C~  tc)  3 i ~ *£>J  + 

where  the  suffixes  1 and  2 indicate  the  direction  of  the  increments. 

The  computation  of  the  increments  requires  the  adjustment  of  the  slip  line 
elements  ds^  and  ds2  so  that  Eqs.  (39)  and  (4C)  are  satisfied. 

B.  Graphical  finite  difference 

This  method  was  developed  by  De  Josselin  de  Jong  £-3^  and  applied  to 
problems  of  failure  in  a two-dimensional  medium  in  plane  strain. 

r 

Use  is  made  of  the  properties  of  the  pole  trail  in  the  stress  plane;  yielding 
an  additional  set  of  equations .. 

The  coordinates  of  the  pole  P are: 

<%■=.  //+ 
z 3 


m 


Eqs.  (36)  and  (37)  require  multiplication  by  r: 

SinX.^ 

The  partial  derivatives  of  (t.££)  and  (t^)  with  respect  to  sA  (i  = 1,  2) 
2£_£1  a //-*.  (41) 

cYt-xa)  «.  (42) 

Multiply  (41)  by  cos  2^  and  (42)  by  sin2^and  add: 


l/rezy*' '***/' 


(43) 


Multiply  (41)  by  (-sin^.sinZyO  and  (42)  by  (1  + siny?.cos2^)  and  add: 

1 . X-3T>  ^ 


we  obtain: 


By  substitution  oi  the  right  hand  terms  of: 

» i'+l**' -?/*•) 

- / -^r'-YA) 


^sv*» y?  +-  toii'j' J * 4 ^°yv/ • 
Multiplication  of  Eq.  (43)  by  cosy*  yields: 


*^<h)  * ^<45> 


and  from  Eq.  (44)  we  obtain: 


(46) 


are*. 


By  adding  Eqs,  (45)  and  (46)  we  obtain  the  left  hand  term  of  Eq.  (36);  by  subtraction, 
the  left  hand  expression  of  Eq.  (37) . Then  introducing  for  (2^-  + y> ) = (2/C?2  + 90°) 
and  tor  (2  f —f)  = (2/3  ^ -90°) , we  obtain: 

'‘"fi ‘ST /*'*>)  “ . rri.eoz/3,)  (48) 

Dividing  through  by  cos/3  and  cos/3  respectively,  and  rearranging: 


^ £-*W- 


In  finite  difference  form: 


2^  = - Zfd*' . co^Jj.(irz).(^sr  ie*/*,)  (51) 

= 7W, ^(r.^X - ^ry  (52) 

The  pole  trail  method  combines  the  relationships  between  the  quantities 
(r.v^,  ) and  y<,  and  the  lengths  of  the  slip  line  elements  ds^  and  ds2  in  their 

repsective  directions^  and  accor<^n9  Eqs*  (56)  an<^  (57) ' 
quantities  (r .T^)  and  (r.C^  ) with  respect  to  the  same  elements,  according  to 
Eqs.  (49)  and  (50). 

This  combination  of  the, geometrical  and  physical  relationships  can  be 
used  to  advantage  in  two-dimensional , plane  strain  problems  to  construct  the 

•lip  line  fields  geometrically  by  finite  differences. 

The  method  was  applied  to  the  problem  under  consideration  but  failed  to 
give  satisfactory  solutions,  because  of  the  arbitrariness  of  the  choice  between  the 

49398 


i 
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terms  and  m coys  on  the  right  hand  side  of  Eqs.  (36)  and  (37)  and  of 

Eqs.  (49)  and  (50). 

The  attempt  showed,  that  while  the  grahpical  solution  is  theoretically  valid, 
it  is  impractical  for  solutions  involving  radial  distributions  of  stresses  in  three- 
dimensional  problems. 

A numerical  finite  difference  method  was  therefore  indicated. 

C.  Numerical  finite  difference  method 

The  differential  equations  of  limit  equilibrium  were  given  in  (36)  and  (37). 
Assuming  that  point  1 and  point  2 in  Fig.  11  are  points  of  the  slip  line 
field  with  known  quantities: 

r, , <£,  * % > 

and:  V ^ "A  * ^ 

Eq.  (36)  applies  to  point  1 located  on  the  Sj-slip  line. 

Eq.  (37)  applies  to  point  2 located  on  the  s2-slip  line. 

We  provide  the  relevant  quantities  with  a second  suffix  to  indicate  the 

location,  of  the  points. 

At  point  1 (Eq.  36): 


eoi'f-\r-  A S)*1  A S) 

At  point  2 (Eq.  (37): 


s - A)* 0/ "/V  ‘ ^ A,*) 


(53) 


y l_  A.s.  (M) 


Point  3 is  reached  along  the  s^slip  line  through  point  1 and  along  the  saline 
through  point  2.  Omitting  the  second  suffix  of  2s  , because  ^ is  related  to 


a 


point  1 and  3.5^  to  point  2,  Eqs.  (53)  and  (54)  are  in  finite  difference  form: 


<°W-  ‘i/vSl-i 

' eisx  y « s'i 


Since: 


^ (*r v*,,)  - *j ■ -z/.C^>/ 

d.  ~ 'Zj  ■ v/nj  - 

d+<  •<&-<& 

<*¥k  * <%-  t< 

Eqs.  (55)  and  (56)  can  be  written  in  the  following  form: 


(55) 


(56) 


0Qif{^j  * -fas*'’  \,+ 


f(h  s - v s)-  //&  -A)  -+/%■*“ n%*  ^ ^ 


tOi 


(57) 

(58) 


The  unknowns  in  these  two  equations  are  VS  ' ^ ' dsiandds2* 

Three  additional  equations  are  needed  to  solve  for  the  unknowns . In  order  to  ob- 
tain these  additional  relationships  the  slip  lines  will  be  considered  as  a sequence  of 
circular  arc  elements  of  varying  curvature. 


When  the  radius  of  curvature  at  a point  of  the  s^-slip  line  is  R^»  the  in 

p 

crease  of  ip  in  the  direction  of  the  slip  line,  d^,  corresponds  to  a slip  line 
distance,  dsA  = Rjd^.  An  identical  relationship  applies  to  the  s2~slip  line; 

“*2  * V*/v 


This  property  will  be  used  to  establish  the  coordinates  of  the  intersection 
x.r.t  3 of  the  s -slip  line  and  s„-slip  line  through  the  points  i.  and  2.  respec- 


Fig.  11  shows  how  these  coordinates  can  be  computed: 


+ X ds>  (***(*'.'  r 


*i  = ^ £ A/'*'  "/V  + (62) 

From:  ^ = yt'-h  f - 'j'r^ 

the  trigonometric  expressions  in/3  can  be  replaced  by  expressions  in  */'' 

A-  - * - ^ •>/*<.,•  *** 

A*-  * ;/A>  * ^ 

3y  substitutions  in  Eqs.  (57)  to  (62)  inclusive  we  obtain  the  six  equations: 


‘yVvS-V^,)WvS ^tM) 


t, « \ +• 





1 


= */  t + rfj/ 


with  the  six  unknown  quantities:  , ds^  and  ds2 . 

We  can  therefore  solve  for  the  unknowns.  The  inclusion,  however,  of  the 
trigonometric  functions  of  requires  a solution  by  numerical  computations. 

The  choice  of  the  intermediate  principal  stress  G~t 


In  Chapter  II.  A we  discussed  the  indeterminedness  of  the  tangential  stress 
(f  . In  Eq.  (28)  the  factor  oC  expresses  this  fact. 

The  value  of  this  factor  can  be  determined  only  at  the  axis  of  symmetry 
and  at  the  boundary  of  the  cavity,  where  the  quantities  r,  G^,  y , /3^,  and 
are  known. 

Along  the  axis  of  symmetry  the  following  quantities  are  known: 

Z s O 

ft  » f *+(**-£} 

/»*  » l**-  t) 

Rewriting  Eqs.  (36)  and  (37)  in  terms  of  the  directional  derivatives: 

+ x?gy>.(ir.)  ^3gi 


J 


Introduction  of  the  known  quantities  of  r,  y- , /3^  and  /3  , yields: 

*0*  y>.  C£.  jg  = _ + W . S*  n ^ . s i V* 

OJ  y?.  - — -f-  4.  oC . , l*ir> 

*~"V*  O 

Dividing  through  by  '7'  and  putting  = cos  /3,  and  = cos  / ?0 

m ^ / 1 *>s„  / ^ 


yields: 


£OJ>  y/.  Ceifi'  at  — 4-  d.3  i-r>  y)  ' * * >V3.i 

cftf-s  yl.  io*/3z  i 4-  s<y'fQi 

Since:  ^ •*  ^oTS’J 

A - /A-  0 

and:  y » (Is  - ‘-‘) 

we  obtain  one  equation: 

4o-styo-i.S‘J.*‘*  «T_  - o4.  iS^. 


from  which: 


hence; 


c?<  - _ / 


/ ,ir  I'uj 

Ut 

The  condition  of  failure  and  the  equilibrium  conditions  along  the  axis  of 
symmetry  require  an  intermediate  principal  stress  of  the  same  magnitude  as  the 

minor  principal  stress  along  that  aixs. 

The  numerical  computations  of  the  failure  quantities  in  the  stress  field 
were  carried  out  on  the  assumption  that  ci  = -1  throughout  the  medium. 

Experimental  evidence  is  required  to  prove  the  validity  of  this  assumption. 
Numerical  computations  with  values  of  -1<*<+1  will  *how  «•»  significance  of  varl- 
ations  in  the  intermediate  principal  stress. 


J/ 


Chapter  III 

Conclusions  concerning  the  feasibility  of  mass  earth  movement  by 

cavity  pressure. 

The  investigation  has  shown  the  feasibility  of  mass  earth  movement  by 
.r.itiating  a failure  condition  through  a rapid  build-up  of  liquid  pressure  inside  a 
cav.ty  at  some  distance  from  the  free  boundary  of  a soil  medium. 

The  magnitude  of  the  critical  pressure  level  depends  primarily  on  the  dis- 
tance of  the  centre  of  pressure  from  the  boundary  and  on  the  mechanical  properties 
cf  the  soil.  The  larger  the  cavity,  the  smaller  the  pressure  required  to  cause 
failure . 

In  a dry  granular  material  (dry  sand)  water  will  penetrate  into  the  soil  by 
the  pressure  gradient  and  saturate  the  immediate  vicinity  of  the  cavity.  The 
saturated  shell  constitutes  the  intermediary  for  the  transfer  of  the  cavity  pressure 
to  the  surrounding  soil. 

In  a saturated  granular  material  a rapid  buiid-up  of  the  liquid  pressure 
inside  the  cavity  propagates  through  the  pore-water  and  reaches  the  watertable 
without  noticeable  delay.  As  a consequence  a critical  flow  gradient  inside  the 
soil  medium  develops  and  while  the  entire  soil  mass  becomes  involved  in  a 
failure  condition,  erosion  of  the  cohesionless  material  is  initiated  at  the  boundary. 

In  a saturated  or  dry  cohesive  material  the  propagation  of  the  cavity 

t 

pressure  depends  on  the  nature  of  the  clay  material.  When  the  void  channel 
system  constitutes  a "closed"  .system,  the  cavity  pressure  will  be  transferred 
mostly  to  the  immediate  vicinity  of  the  cavity  and  failure  will  occur  as  if  the 
material  would  be  an  impervious  medium. 


in  materials  with  an  "open"  void  channel  system  the  cavity  pressure  creates 
.condition  comparable  to  that  of  a saturated  granular  medium. 

In  all  soil  media  the  state  of  failure  by  rapid  pressure  build-up  is  followed 
deterioration  of  the  material.  The  time  period  for  maintaining  critically  high 
pressures  has  shown  to  be  comparitively  shorten  the  order  of  seconds).  The 
vc. -me  displacements  as  a result  of  the  expansion  of  the  cavity,  preceding  the 
state  of  failure,  are  compantive*y  small.  In  the  first  phase  of  the  treatment  high 
pressure  input  and  a small  water  quantity  are  required.  In  the  second  phase  the 
pressure  cun  be  decreased,  but  a larger  quantity  of  water  is  required  to  maintain 
c liquid  flow  of  soil  particles  and  deteriorated  soil  materials. 

» 

Example  of  a continuously  operating  mass  earth  movement  system 

As  an  example  of  a continuously  operating  system.  Fig.  11  shows  a diagram 
of  the  mechanism  required  to  move  the  soil  in  excavating  a canal  profile. 

A liquid  pressure  feeding  line  is  located  across  the  canal  axis  and  connected 
with  a number  of  casings,  which  are  gradually  lowered  into  the  soil,  wnile 
pressure  is  constantly  applied.  The  liquid  flow  moves  the  soil  out  and  a system 
01  suction  lines  at  the  toe  line  removes  the  suspended  material. 

This  material  can  be  pumped  over  the  sides  and  deposited  alongside  at 
any  suitable  location. 

The  excavation  is  supposed  to  remain  in  open  connection  with  a reservoir 
(river,  ocean)  allowing  for  a constant,  unlimited  supply  of  water. 

Water  is  used  both  as  a pressure  medium  and  a transportation  device. 

*he  entire  mechanism  can  be  moved  along  the  canal  trace,  while  the  liquid 


material  is  deposited  alongside. 


By  adding  suitable  chemical  admixtures  the  deterioration  of  cohesive 
material  can  be  speeded  up.  When  pseudo-rock,  materials  are  encountered  the 
initiation  of  a failure  condition,  using  liquid  pressure,  can  be  promoted  by 
..sing  suitable  explosives  inside  the  cavity.  In  granular  materials  and  cohesive 
soils  of  medium  strength  such  measures  will  not  be  necessary,  provided  that 
an  average  depth  of  the  cavity  is  maintained. 

In  cohesionless,  dry  granular  materials  the  addition  of  water  to  the 
soil  promotes  the  occurrence  of  flow  slides,  which-when  kept  in  continuous 
motion  by  transportation  through  the  suction  lines-do  require  little  energy. 

In  cohesionless,  wet  granular  material  the  removal  of  material  along 
the  toe  of  the  slope  will  generally  be  sufficient  to  cause  flow  slides  as  a con- 
tinuous earth  moving  mechanism. 
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